INHBA is Enriched in HPV-negative Oropharyngeal Squamous Cell Carcinoma and Promotes Cancer Progression.

Patients with oropharyngeal squamous cell carcinoma (OPSCC) caused by human papilloma virus (HPV) exhibit a better prognosis than those with HPV-negative OPSCC. This study investigated the distinct molecular pathways that delineate HPV-negative from HPV-positive OPSCC to identify biologically relevant therapeutic targets. Bulk mRNA from 23 HPV-negative and 39 HPV-positive OPSCC tumors (n = 62) was sequenced to uncover the transcriptomic profiles. Differential expression followed by gene set enrichment analysis was performed to outline the top enriched biological process in the HPV-negative compared with HPV-positive entity. INHBA, the highest overexpressed gene in the HPV-negative tumor, was knocked down. Functional assays (migration, proliferation, cell death, stemness) were conducted to confirm the target's oncogenic role. Correlation analyses to reveal its impact on the tumor microenvironment were performed. We revealed that epithelial-to-mesenchymal transition (EMT) is the most enriched process in HPV-negative compared with HPV-positive OPSCC, with INHBA (inhibin beta A subunit) being the top upregulated gene. INHBA knockdown downregulated the expression of EMT transcription factors and attenuated migration, proliferation, stemness, and cell death resistance of OPSCC cells. We uncovered that INHBA associates with a pro-tumor microenvironment by negatively correlating with antitumor CD8+ T and B cells while positively correlating with pro-tumor M1 macrophages. We identified three miRNAs that are putatively involved in repressing INHBA expression. Our results indicate that the upregulation of INHBA is tumor-promoting. We propose INHBA as an attractive therapeutic target for the treatment of INHBA-enriched tumors in patients with HPV-negative OPSCC to ameliorate prognosis. SIGNIFICANCE: Patients with HPV-negative OPSCC have a poorer prognosis due to distinct molecular pathways. This study reveals significant transcriptomic differences between HPV-negative and HPV-positive OPSCC, identifying INHBA as a key upregulated gene in HPV-negative OPSCC's oncogenic pathways. INHBA is crucial in promoting EMT, cell proliferation, and an immunosuppressive tumor environment, suggesting its potential as a therapeutic target for HPV-negative OPSCC.

Acetylation-induced proteasomal degradation of the activated glucocorticoid receptor limits hormonal signaling.

Glucocorticoid (GC) signaling is essential for mounting a stress response, however, chronic stress or prolonged GC therapy downregulates the GC receptor (GR), leading to GC resistance. Regulatory mechanisms that refine this equilibrium are not well understood. Here, we identify seven lysine acetylation sites in the amino terminal domain of GR, with lysine 154 (Lys154) in the AF-1 region being the dominant acetyl-acceptor. GR-Lys154 acetylation is mediated by p300/CBP in the nucleus in an agonist-dependent manner and correlates with transcriptional activity. Deacetylation by NAD+-dependent SIRT1 facilitates dynamic regulation of this mark. Notably, agonist-binding to both wild-type GR and an acetylation-deficient mutant elicits similar short-term target gene expression. In contrast, upon extended treatment, the polyubiquitination of the acetylation-deficient GR mutant is impaired resulting in higher protein stability, increased chromatin association and prolonged transactivation. Taken together, reversible acetylation fine-tunes duration of the GC response by regulating proteasomal degradation of activated GR.

Protecting RNA quality for spatial transcriptomics while improving immunofluorescent staining quality.

In comparison to bulk sequencing or single cell sequencing, spatial transcriptomics preserves the spatial information in tissue slices and can even be mapped to immunofluorescent stainings, allowing translation of gene expression information into their spatial context. This enables to unravel complex interactions of neighboring cells or to link cell morphology to transcriptome data. The 10× Genomics Visium platform offers to combine spatial transcriptomics with immunofluorescent staining of cryo-sectioned tissue slices. We applied this technique to fresh frozen mouse brain slices and developed a protocol that still protects RNA quality while improving buffers for immunofluorescent staining. We investigated the impact of various parameters, including fixation time and buffer composition, on RNA quality and antibody binding. Here, we propose an improved version of the manufacturer protocol, which does not alter RNA quality and facilitates the use of multiple additional antibodies that were not compatible with the manufacturer protocol before. Finally, we discuss the influence of various staining parameters, which contribute to the development of application specific staining protocols.

Gene regulatory network inference with popInfer reveals dynamic regulation of hematopoietic stem cell quiescence upon diet restriction and aging.

Inference of gene regulatory networks (GRNs) can reveal cell state transitions from single-cell genomics data. However, obstacles to temporal inference from snapshot data are difficult to overcome. Single-nuclei multiomics data offer means to bridge this gap and derive temporal information from snapshot data using joint measurements of gene expression and chromatin accessibility in the same single cells. We developed popInfer to infer networks that characterize lineage-specific dynamic cell state transitions from joint gene expression and chromatin accessibility data. Benchmarking against alternative methods for GRN inference, we showed that popInfer achieves higher accuracy in the GRNs inferred. popInfer was applied to study single-cell multiomics data characterizing hematopoietic stem cells (HSCs) and the transition from HSC to a multipotent progenitor cell state during murine hematopoiesis across age and dietary conditions. From networks predicted by popInfer, we discovered gene interactions controlling entry to/exit from HSC quiescence that are perturbed in response to diet or aging.

The HLA ligandome of oropharyngeal squamous cell carcinomas reveals shared tumour-exclusive peptides for semi-personalised vaccination.

BACKGROUND: The immune peptidome of OPSCC has not previously been studied. Cancer-antigen specific vaccination may improve clinical outcome and efficacy of immune checkpoint inhibitors such as PD1/PD-L1 antibodies. METHODS: Mapping of the OPSCC HLA ligandome was performed by mass spectrometry (MS) based analysis of naturally presented HLA ligands isolated from tumour tissue samples (n = 40) using immunoaffinity purification. The cohort included 22 HPV-positive (primarily HPV-16) and 18 HPV-negative samples. A benign reference dataset comprised of the HLA ligandomes of benign haematological and tissue datasets was used to identify tumour-associated antigens. RESULTS: MS analysis led to the identification of naturally HLA-presented peptides in OPSCC tumour tissue. In total, 22,769 peptides from 9485 source proteins were detected on HLA class I. For HLA class II, 15,203 peptides from 4634 source proteins were discovered. By comparative profiling against the benign HLA ligandomic datasets, 29 OPSCC-associated HLA class I ligands covering 11 different HLA allotypes and nine HLA class II ligands were selected to create a peptide warehouse. CONCLUSION: Tumour-associated peptides are HLA-presented on the cell surfaces of OPSCCs. The established warehouse of OPSCC-associated peptides can be used for downstream immunogenicity testing and peptide-based immunotherapy in (semi)personalised strategies.

Increased longevity due to sexual activity in mole-rats is associated with transcriptional changes in the HPA stress axis.

Sexual activity and/or reproduction are associated with a doubling of life expectancy in the long-lived rodent genus Fukomys. To investigate the molecular mechanisms underlying this phenomenon, we analyzed 636 RNA-seq samples across 15 tissues. This analysis suggests that changes in the regulation of the hypothalamic-pituitary-adrenal stress axis play a key role regarding the extended life expectancy of reproductive vs. non-reproductive mole-rats. This is substantiated by a corpus of independent evidence. In accordance with previous studies, the up-regulation of the proteasome and so-called 'anti-aging molecules', for example, dehydroepiandrosterone, is linked with enhanced lifespan. On the other hand, several of our results are not consistent with knowledge about aging of short-lived model organisms. For example, we found the up-regulation of the insulin-like growth factor 1/growth hormone axis and several other anabolic processes to be compatible with a considerable lifespan prolongation. These contradictions question the extent to which findings from short-lived species can be transferred to longer-lived ones.

Correction: The acetyltransferase GCN5 maintains ATRA-resistance in non-APL AML.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The acetyltransferase GCN5 maintains ATRA-resistance in non-APL AML.

To date, only one subtype of acute myeloid leukemia (AML), acute promyelocytic leukemia (APL) can be effectively treated by differentiation therapy utilizing all-trans retinoic acid (ATRA). Non-APL AMLs are resistant to ATRA. Here we demonstrate that the acetyltransferase GCN5 contributes to ATRA resistance in non-APL AML via aberrant acetylation of histone 3 lysine 9 (H3K9ac) residues maintaining the expression of stemness and leukemia associated genes. We show that inhibition of GCN5 unlocks an ATRA-driven therapeutic response. This response is potentiated by coinhibition of the lysine demethylase LSD1, leading to differentiation in most non-APL AML. Induction of differentiation was not correlated to a specific AML subtype, cytogenetic, or mutational status. Our study shows a previously uncharacterized role of GCN5 in maintaining the immature state of leukemic blasts and identifies GCN5 as a therapeutic target in AML. The high efficacy of the combined epigenetic treatment with GCN5 and LSD1 inhibitors may enable the use of ATRA for differentiation therapy of non-APL AML. Furthermore, it supports a strategy of combined targeting of epigenetic factors to improve treatment, a concept potentially applicable for a broad range of malignancies.

Analysis of the coding sequences of clownfish reveals molecular convergence in the evolution of lifespan.

BACKGROUND: Standard evolutionary theories of aging postulate that reduced extrinsic mortality leads to evolution of longevity. Clownfishes of the genus Amphiprion live in a symbiotic relationship with sea anemones that provide protection from predators. We performed a survey and identified at least two species with a lifespan of over 20 years. Given their small size and ease of captive reproduction, clownfish lend themselves as experimental models of exceptional longevity. To identify genetic correlates of exceptional longevity, we sequenced the transcriptomes of Amphiprion percula and A. clarkii and performed a scan for positively-selected genes (PSGs). RESULTS: The PSGs that we identified in the last common clownfish ancestor were compared with PSGs detected in long-lived mole rats and short-lived killifishes revealing convergent evolution in processes such as mitochondrial biogenesis. Among individual genes, the Mitochondrial Transcription Termination Factor 1 (MTERF1), was positively-selected in all three clades, whereas the Glutathione S-Transferase Kappa 1 (GSTK1) was under positive selection in two independent clades. For the latter, homology modelling strongly suggested that positive selection targeted enzymatically important residues. CONCLUSIONS: These results indicate that specific pathways were recruited in independent lineages evolving an exceptionally extended or shortened lifespan and point to mito-nuclear balance as a key factor.

Higher gene expression stability during aging in long-lived giant mole-rats than in short-lived rats.

Many aging-associated physiological changes are known to occur in short- and long-lived species with different trajectories. Emerging evidence suggests that numerous life history trait differences between species are based on interspecies variations in gene expression. Little information is available, however, about differences in transcriptome changes during aging between mammals with diverging lifespans. For this reason, we studied the transcriptomes of five tissue types and two age cohorts of two similarly sized rodent species with very different lifespans: laboratory rats (Rattus norvegicus) and giant mole-rats (Fukomys mechowii), with maximum lifespans of 3.8 and more than 20 years, respectively. Our findings show that giant mole-rats exhibit higher gene expression stability during aging than rats. Although well-known aging signatures were detected in all tissue types of rats, they were found in only one tissue type of giant mole-rats. Furthermore, many differentially expressed genes that were found in both species were regulated in opposite directions during aging. This suggests that expression changes which cause aging in short-lived species are counteracted in long-lived species. Taken together, we conclude that expression stability in giant mole rats (and potentially in African mole-rats in general) may be one key factor for their long and healthy life.

Species comparison of liver proteomes reveals links to naked mole-rat longevity and human aging.

BACKGROUND: Mammals display a wide range of variation in their lifespan. Investigating the molecular networks that distinguish long- from short-lived species has proven useful to identify determinants of longevity. Here, we compared the livers of young and old long-lived naked mole-rats (NMRs) and the phylogenetically closely related, shorter-lived, guinea pigs using an integrated omics approach. RESULTS: We found that NMR livers display a unique expression pattern of mitochondrial proteins that results in distinct metabolic features of their mitochondria. For instance, we observed a generally reduced respiration rate associated with lower protein levels of respiratory chain components, particularly complex I, and increased capacity to utilize fatty acids. Interestingly, we show that the same molecular networks are affected during aging in both NMRs and humans, supporting a direct link to the extraordinary longevity of both species. Finally, we identified a novel detoxification pathway linked to longevity and validated it experimentally in the nematode Caenorhabditis elegans. CONCLUSIONS: Our work demonstrates the benefits of integrating proteomic and transcriptomic data to perform cross-species comparisons of longevity-associated networks. Using a multispecies approach, we show at the molecular level that livers of NMRs display progressive age-dependent changes that recapitulate typical signatures of aging despite the negligible senescence and extraordinary longevity of these rodents.

Naked mole-rat transcriptome signatures of socially suppressed sexual maturation and links of reproduction to aging.

BACKGROUND: Naked mole-rats (NMRs) are eusocially organized in colonies. Although breeders carry the additional metabolic load of reproduction, they are extremely long-lived and remain fertile throughout their lifespan. This phenomenon contrasts the disposable soma theory of aging stating that organisms can invest their resources either in somatic maintenance, enabling a longer lifespan, or in reproduction, at the cost of longevity. Here, we present a comparative transcriptome analysis of breeders vs. non-breeders of the eusocial, long-lived NMR vs. the polygynous and shorter-lived guinea pig (GP). RESULTS: Comparative transcriptome analysis of tissue samples from ten organs showed, in contrast to GPs, low levels of differentiation between sexes in adult NMR non-breeders. After transition into breeders, NMR transcriptomes are markedly sex-specific, show pronounced feedback signaling via gonadal steroids, and have similarities to reproductive phenotypes in African cichlid fish, which also exhibit social status changes between dominant and subordinate phenotypes. Further, NMRs show functional enrichment of status-related expression differences associated with aging. Lipid metabolism and oxidative phosphorylation-molecular networks known to be linked to aging-were identified among most affected gene sets. Remarkably and in contrast to GPs, transcriptome patterns associated with longevity are reinforced in NMR breeders. CONCLUSION: Our results provide comprehensive and unbiased molecular insights into interspecies differences between NMRs and GPs, both in sexual maturation and in the impact of reproduction on longevity. We present molecular evidence that sexual maturation in NMRs is socially suppressed. In agreement with evolutionary theories of aging in eusocial organisms, we have identified transcriptome patterns in NMR breeders that-in contrast to the disposable soma theory of aging-may slow down aging rates and potentially contribute to their exceptional long life- and healthspan.

Long-lived rodents reveal signatures of positive selection in genes associated with lifespan.

The genetics of lifespan determination is poorly understood. Most research has been done on short-lived animals and it is unclear if these insights can be transferred to long-lived mammals like humans. Some African mole-rats (Bathyergidae) have life expectancies that are multiple times higher than similar sized and phylogenetically closely related rodents. To gain new insights into genetic mechanisms determining mammalian lifespans, we obtained genomic and transcriptomic data from 17 rodent species and scanned eleven evolutionary branches associated with the evolution of enhanced longevity for positively selected genes (PSGs). Indicating relevance for aging, the set of 250 identified PSGs showed in liver of long-lived naked mole-rats and short-lived rats an expression pattern that fits the antagonistic pleiotropy theory of aging. Moreover, we found the PSGs to be enriched for genes known to be related to aging. Among these enrichments were "cellular respiration" and "metal ion homeostasis", as well as functional terms associated with processes regulated by the mTOR pathway: translation, autophagy and inflammation. Remarkably, among PSGs are RHEB, a regulator of mTOR, and IGF1, both central components of aging-relevant pathways, as well as genes yet unknown to be aging-associated but representing convincing functional candidates, e.g. RHEBL1, AMHR2, PSMG1 and AGER. Exemplary protein homology modeling suggests functional consequences for amino acid changes under positive selection. Therefore, we conclude that our results provide a meaningful resource for follow-up studies to mechanistically link identified genes and amino acids under positive selection to aging and lifespan determination.

Tissue-, sex-, and age-specific DNA methylation of rat glucocorticoid receptor gene promoter and insulin-like growth factor 2 imprinting control region.

Tissue-, sex-, and age-specific epigenetic modifications such as DNA methylation are largely unknown. Changes in DNA methylation of the glucocorticoid receptor gene (NR3C1) and imprinting control region (ICR) of IGF2 and H19 genes during the lifespan are particularly interesting since these genes are susceptible to epigenetic modifications by prenatal stress or malnutrition. They are important regulators of development and aging. Methylation changes of NR3C1 affect glucocorticoid receptor expression, which is associated with stress sensitivity and stress-related diseases predominantly occurring during aging. Methylation changes of IGF2/H19 affect growth trajectory and nutrient use with risk of metabolic syndrome. Using a locus-specific approach, we characterized DNA methylation patterns of different Nr3c1 promoters and Igf2/H19 ICR in seven tissues of rats at 3, 9, and 24 mo of age. We found a complex pattern of locus-, tissue-, sex-, and age-specific DNA methylation. Tissue-specific methylation was most prominent at the shores of the Nr3c1 CpG island (CGI). Sex-specific differences in methylation peaked at 9 mo. During aging, Nr3c1 predominantly displayed hypomethylation mainly in females and at shores, whereas hypermethylation occurred within the CGI. Igf2/H19 ICR exhibited age-related hypomethylation occurring mainly in males. Methylation patterns of Nr3c1 in the skin correlated with those in the cortex, hippocampus, and hypothalamus. Skin may serve as proxy for methylation changes in central parts of the hypothalamic-pituitary-adrenal axis and hence for vulnerability to stress- and age-associated diseases. Thus, we provide in-depth insight into the complex DNA methylation changes of rat Nr3c1 and Igf2/H19 during aging that are tissue and sex specific.

PosiGene: automated and easy-to-use pipeline for genome-wide detection of positively selected genes.

Many comparative genomics studies aim to find the genetic basis of species-specific phenotypic traits. A prevailing strategy is to search genome-wide for genes that evolved under positive selection based on the non-synonymous to synonymous substitution ratio. However, incongruent results largely due to high false positive rates indicate the need for standardization of quality criteria and software tools. Main challenges are the ortholog and isoform assignment, the high sensitivity of the statistical models to alignment errors and the imperative to parallelize large parts of the software. We developed the software tool PosiGene that (i) detects positively selected genes (PSGs) on genome-scale, (ii) allows analysis of specific evolutionary branches, (iii) can be used in arbitrary species contexts and (iv) offers visualization of the results for further manual validation and biological interpretation. We exemplify PosiGene's performance using simulated and real data. In the simulated data approach, we determined a false positive rate <1%. With real data, we found that 68.4% of the PSGs detected by PosiGene, were shared by at least one previous study that used the same set of species. PosiGene is a user-friendly, reliable tool for reproducible genome-wide identification of PSGs and freely available at https://github.com/gengit/PosiGene.

Parallel evolution of genes controlling mitonuclear balance in short-lived annual fishes.

The current molecular understanding of the aging process derives almost exclusively from the study of random or targeted single-gene mutations in highly inbred laboratory species, mostly invertebrates. Little information is available as to the genetic mechanisms responsible for natural lifespan variation and the evolution of lifespan, especially in vertebrates. Here, we investigated the pattern of positive selection in annual (i.e., short-lived) and nonannual (i.e., longer-lived) African killifishes to identify a genomic substrate for evolution of annual life history (and reduced lifespan). We identified genes under positive selection in all steps of mitochondrial biogenesis: mitochondrial (mt) DNA replication, transcription from mt promoters, processing and stabilization of mt RNAs, mt translation, assembly of respiratory chain complexes, and electron transport chain. Signs of paralleled evolution (i.e., evolution in more than one branch of Nothobranchius phylogeny) are observed in four out of five steps. Moreover, some genes under positive selection in Nothobranchius are under positive selection also in long-lived mammals such as bats and mole-rats. Complexes of the respiratory chain are formed in a coordinates multistep process where nuclearly and mitochondrially encoded components are assembled and inserted into the inner mitochondrial membrane. The coordination of this process is named mitonuclear balance, and experimental manipulations of mitonuclear balance can increase longevity of laboratory species. Our data strongly indicate that these genes are also casually linked to evolution lifespan in vertebrates.

FRAMA: from RNA-seq data to annotated mRNA assemblies.

BACKGROUND: Advances in second-generation sequencing of RNA made a near-complete characterization of transcriptomes affordable. However, the reconstruction of full-length mRNAs via de novo RNA-seq assembly is still difficult due to the complexity of eukaryote transcriptomes with highly similar paralogs and multiple alternative splice variants. Here, we present FRAMA, a genome-independent annotation tool for de novo mRNA assemblies that addresses several post-assembly tasks, such as reduction of contig redundancy, ortholog assignment, correction of misassembled transcripts, scaffolding of fragmented transcripts and coding sequence identification. RESULTS: We applied FRAMA to assemble and annotate the transcriptome of the naked mole-rat and assess the quality of the obtained compilation of transcripts with the aid of publicy available naked mole-rat gene annotations. Based on a de novo transcriptome assembly (Trinity), FRAMA annotated 21,984 naked mole-rat mRNAs (12,100 full-length CDSs), corresponding to 16,887 genes. The scaffolding of 3488 genes increased the median sequence information 1.27-fold. In total, FRAMA detected and corrected 4774 misassembled genes, which were predominantly caused by fusion of genes. A comparison with three different sources of naked mole-rat transcripts reveals that FRAMA's gene models are better supported by RNA-seq data than any other transcript set. Further, our results demonstrate the competitiveness of FRAMA to state of the art genome-based transcript reconstruction approaches. CONCLUSION: FRAMA realizes the de novo construction of a low-redundant transcript catalog for eukaryotes, including the extension and refinement of transcripts. Thereby, results delivered by FRAMA provide the basis for comprehensive downstream analyses like gene expression studies or comparative transcriptomics. FRAMA is available at https://github.com/gengit/FRAMA .

Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish.

The killifish Nothobranchius furzeri is the shortest-lived vertebrate that can be bred in the laboratory. Its rapid growth, early sexual maturation, fast aging, and arrested embryonic development (diapause) make it an attractive model organism in biomedical research. Here, we report a draft sequence of its genome that allowed us to uncover an intra-species Y chromosome polymorphism representing-in real time-different stages of sex chromosome formation that display features of early mammalian XY evolution "in action." Our data suggest that gdf6Y, encoding a TGF-ß family growth factor, is the master sex-determining gene in N. furzeri. Moreover, we observed genomic clustering of aging-related genes, identified genes under positive selection, and revealed significant similarities of gene expression profiles between diapause and aging, particularly for genes controlling cell cycle and translation. The annotated genome sequence is provided as an online resource (http://www.nothobranchius.info/NFINgb).

Unusual ratio between free thyroxine and free triiodothyronine in a long-lived mole-rat species with bimodal ageing.

Ansell's mole-rats (Fukomys anselli) are subterranean, long-lived rodents, which live in eusocial families, where the maximum lifespan of breeders is twice as long as that of non-breeders. Their metabolic rate is significantly lower than expected based on allometry, and their retinae show a high density of S-cone opsins. Both features may indicate naturally low thyroid hormone levels. In the present study, we sequenced several major components of the thyroid hormone pathways and analyzed free and total thyroxine and triiodothyronine in serum samples of breeding and non-breeding F. anselli to examine whether a) their thyroid hormone system shows any peculiarities on the genetic level, b) these animals have lower hormone levels compared to euthyroid rodents (rats and guinea pigs), and c) reproductive status, lifespan and free hormone levels are correlated. Genetic analyses confirmed that Ansell's mole-rats have a conserved thyroid hormone system as known from other mammalian species. Interspecific comparisons revealed that free thyroxine levels of F. anselli were about ten times lower than of guinea pigs and rats, whereas the free triiodothyronine levels, the main biologically active form, did not differ significantly amongst species. The resulting fT4:fT3 ratio is unusual for a mammal and potentially represents a case of natural hypothyroxinemia. Comparisons with total thyroxine levels suggest that mole-rats seem to possess two distinct mechanisms that work hand in hand to downregulate fT4 levels reliably. We could not find any correlation between free hormone levels and reproductive status, gender or weight. Free thyroxine may slightly increase with age, based on sub-significant evidence. Hence, thyroid hormones do not seem to explain the different ageing rates of breeders and non-breeders. Further research is required to investigate the regulatory mechanisms responsible for the unusual proportion of free thyroxine and free triiodothyronine.